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Abstract 

Effects of the 0 parameter are studied in Witten’s model of holographic 4d Yang-Mills, where 
6 is the coefficient of the CP-breaking topological term. First, the gravity background, 
including the full backreaction of the RR form dual to the 6 parameter, is revisited. Then, 
a number of observables are computed holographically: the ground-state energy density, the 
string tension, the ’t Hooft loop, the light scalar glueball mass, the baryon mass scale, the 
critical temperature for deconfinement - and thus the whole (T, 6) phase diagram - and the 
entanglement entropy. A simple rule is provided to derive the 6 corrections to (at least) all 
the CP-neutral observables of the model. Some of the observables we consider can and have 
been in fact studied in pure 4d Yang-Mills on the lattice. In that framework the results, 
obtained in the small 6 regime, are given up to very few powers of 9 2 . The corresponding 
holographic results agree qualitatively with available lattice data and signal an overall mass 
scale reduction by 0. Moreover, being exact in 6 , they provide a benchmark for higher order 
corrections in Yang-Mills. 
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1 Introduction and results 


The top-down prototype of holographic Yang-Mills, due to Witten [Tj, is based on the back¬ 
ground generated by ZM-branes wrapped on a circle with appropriate boundary conditions. 
The held theory dual is a non supersymmetric confining Yang-Mills coupled to massive 
“Kaluza-Klein” (KK) matter in the adjoint representation. In this note we are interested in 
studying the effects of the Yang-Mills 9 angle - the coefficient of the CP-breaking topological 
term in the Lagrangian - in this model. 

There exist interesting studies of the 6 dependence of some observables in lattice Yang- 
Mills, see for example the excellent review [2] and references therein. As usual, these studies 
can be performed in principle at any number of colors N c and for the lattice discretization 
of actual pure Yang-Mills. On the other hand they are challenging, due to a sign problem 
related to the 9 term. For this reason, lattice results are limited to small values of 9 (mainly 
to the first few non trivial orders in the 9 2 expansion around 9 = 0), as they are either 
analytically continued from imaginary 9 values or obtained by computing certain correlators 
at 9 — 0. 

In the holographic model there is no such limitation. When 9 is very small, one can treat 
the corresponding (1-form Ramond-Ramond) field in the dual gravity background as a probe 
[3], It is not difficult, however, to account for the full backreaction of that field and so to 
consider finite generic values of 9. The corresponding gravity solution has been found in 
This is what we will focus on in this note. 

In the holographic framework, the main limitations arise from the fact that the classical 
gravity solution, dual to the field theory, is reliable only in the planar N c ^> 1 limit at “strong 
coupling” A 4 » 1, where A 4 , which can be viewed as the Yang-Mills ’t Hooft coupling at 
the Kaluza-Klein mass scale Mkk , is a parameter which actually measures how much the 
spurious KK matter fields are decoupled from the Yang-Mills theory. When A 4 is large the 
decoupling cannot be achieved. Despite this limit, Witten’s model realizes in a very simple 
way all the crucial expected IR features of the pure Yang-Mills theory. In this note we want 
to provide some novel evidence for this to be true also at finite 0-angle. As we will see, also 
because the 9 term enters in Witten’s Yang-Mills theory in the same way it enters in pure 
Yang-Mills, the holographic approach is able to capture the same qualitative trends of the 
topological effects expected in pure Yang-Mills theory. This justifies the comparison we will 
make with lattice results on the way. 

At large N c and A 4 , the effects of the 9 angle turn out to be actually weighed by the 
combination 0 ~ A 4 9/N c . In order for non-trivial effects due to the 9 angle to be considered, 
the limits have to be taken in such a way that 0 remains finite. 

After having presented the gravity background in Section [2] (see also Appendix [A] for 
a re-derivation), in Section [3] we enroll in the study of the 9 dependence of a number of 
relevant observables in the dual gauge theory. Although the theory is not precisely pure 
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Yang-Mills and it is studied in the planar limit, it allows to derive the exact 9 dependence 
of the observables, providing a complementary view with respect to lattice Yang-Mills. 

Some of the observables we calculate have been studied in lattice Yang-Mills: the fun¬ 
damental string tension, the light scalar glueball mass, and the critical temperature T c for 
deconfinement have all been computed to order o 2 BE]. In all these cases, at the same 
order, the observables in the holographic model agree qualitatively^ with the lattice results 
- namely the 9 2 corrections to the 9 = 0 values are all negative: the 9 term reduces the mass 
scales. 

The ground-state energy density, whose 0(9 2 ) coefficient gives the topological suscepti¬ 
bility, has also been computed on the lattice to order d 4 (the 0(9 6 ) coefficient has been 
considered e.g. in BE]: unfortunately it is not possible to determine its sign yet, as the 
errors are still very large). Also in this case, holographic and lattice data qualitatively agree. 

Interestingly enough, in the holographic model the above d-corrections are just the first 
non trivial terms of (powers of) geometric series in © 2 . More precisely, they come from the 
expansion of functions like (1 + 0 2 ) _a ,a > 0. As such, the holographic model predicts that 
the 0(9 A ) corrections to, say, the string tension and the glueball mass, will be of positive 
sign. Analogously, the 0(9 6 ) coefficient in the ground-state energy density (which has an 
overall negative sign) is predicted to be positive. It would be very interesting to improve 
the lattice precision in order to check whether these predictions are actually realized in pure 
Yang-Mills. 

As another notable feature, the holographic model precisely accounts for the expected 
invariance of the gauge theory observables under 9 —> 9 + 2ti shifts. For instance, the 
expected (see e.g. H) periodic structure of the (T, 9) phase diagram with triple points 
where first-order phase transition lines meet, is explicitly realized. 

In the present work, we also consider the 9 dependence of quantities for which there are 
no lattice results yet, namely the’t Hooft loop (from which we deduce an explicit realization 
of the so-called “oblique confinement”), the mass scale of baryonsj^] which follows the same 
pattern as the observables described above, and the entanglement entropy. The latter is 
studied for two different geometries, the “slab” and the ball. In both cases there is a phase 
transition as the dimension of the entangling region is varied, between a “IR configuration” 
and a “UV configuration” mm- While the former scales like (1 + 0 2 ) “, the latter is 
insensitive to ©, consistently with the fact that the 9 dependence is a IR phenomenon. The 
critical length for the transition grows as (1 + 0 2 ) 1 / 2 . 

All of the © corrections to the observables in the holographic model are extremely simple: 
they are given by multiplying (factors of) the © = 0 result by powers of (1 + © 2 ). It is 

1 'We do not report on a comprehensive quantitative comparison of the holographic model with the lattice 
data. The latter are still affected by large statistical errors in the large N c limit, and the two approaches 
currently seem to present significant discrepancies. 

2 In [TO] baryons have been investigated in the related solution in 0]. 
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actually a-posteriori trivial to guess which power of (1 + © 2 ) appears in each observable. In 
fact, the model has two distinct mass scales: the Kaluza-Klein scale Mkk (which is also the 
glueball mass scale at 0 = 0) and the string tension T s . Their ratio is determined by the 
parameter A 4 ~ T S /M 2 <K . As it turns out, it is sufficient to determine how the two quantities 
Mkk and A 4 (or equivalently T s ) scale with 0 to have the scaling of all the observables. To 
be specific, it is sufficient to include a factor of (1 + @ 2 ) -1 / 2 for each power of Mkk and a 
factor of (1 + 0 2 )” 1 for each power of A 4 appearing in a given observable. 

Although we have no counter-example to this rule, we cannot exclude that observables 
which are sensible to CP parity evade it, due to mixing effects^] Apart from this caveat, 
the rule above would allow to write down the 0 corrections to at least all the CP-neutral 
observables of the model, once their Mkk and A 4 factors are known at 0 = without 
performing any calculation! This is ultimately due to the simple origin of the 6 deformation 
of the background in lid, where it is basically a twist of a two-cycle H31I5]- It would be 
extremely interesting to fold a similar pattern on the lattice, at least for some classes of 
observables. 

As a final aside comment, let us recall that chiral massless fermions, transforming in the 
fundamental of the gauge group, can be introduced in Witten’s model by adding D8-branes, 
giving the prototype “holographic QCD” model of Sakai and Sugimoto [15]. In that setup, 
the axial U(1)a anomaly and the Witten-Veneziano mechanism relating the mass of the 
?/ meson with the topological susceptibility of the unflavored theory are precisely realized. 
At the same time, consistently with what happens in massless QCD, a non zero 6 term in 
the Sakai-Sugimoto model can be rotated away by a chiral rotation of the fermions. As a 
result the topological susceptibility of the flavored theory is zero. In order to see the effects 
of the topological term in this QCD-like setup, one needs to switch to the case where the 
fundamental matter fields are massive. We hope to come back to these issues in the future. 

2 Theta dependence in Holographic Yang-Mills 

The Yang-Mills Euclidean Lagrangian at finite 9 is given b}0 

(2 - i) 

where A = g 2 N c is the ’t Hooft coupling. The CP breaking 6 parameter multiplies the 
topological charge density and as such it behaves like an angle: the physics has to be invariant 

3 The mesonic spectra in m could be of this kind. 

4 This rule would suggest in turn that observables scaling, at 0 = 0, as powers of Aare not 
corrected by 9. 

5 We use the standard normalization conventions Tr [t a tb] = (1/2 )6 a b for the SU(N C ) generators t a . Notice 
moreover that = (1/2 )e^ pa F apa . In Minkowski, the Lagrangian has a minus sign on the first term and 
a plus sign, replacing the imaginary factor —i, on the second term. 
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under shifts 0 —)• 0 + 2n. At the same time, as it is evident from the previous expression, in 
the large N c limit the theory has a non-trivial 0-dependence (as required, e.g. by the large 
N c solution of the U(l )a problem, i.e. why the rf mass is much higher then those of the 
meson octet) only if 9/N c is held fixed. Observables, like the ground-state energy density, 
will thus be functions of 9/N c . 

As it was shown in [TE], a way to reconcile the periodicity requirement and the large 
N c scaling is to allow for the vacuum energy density of the theory to be a multi-branched 
function. Possible vacua, labeled by an integer k, become stable but non-degenerate at 
N c — oc and the ground-state energy, for a given value of 0, is obtained by minimizing, 
with respect to k, a function of (0 + 2nk)/N c . The resulting expression, periodic in 0, is not 
regular everywhere: at certain values of 0 it accounts for a jump between different branches. 

In the holographic Yang-Mills model, this expected behavior is explicitly realized [3]. 

Witten’s original background [T) is sourced by N c 1 ZM-branes (of type IIA string 
theory) wrapped on a circle S X4 , of length 27where periodic (resp. anti-periodic) 
boundary conditions on scalars (resp. fermions) are imposed. In such a way, in the deep IR 
(i.e. at energies E <C M KK ) the original 4 + 1 dimensional theory on the Z04-branes, reduces 
to pure non-supersymmetric SU(N C ) Yang-Mills in 3 + 1 dimensions. In fact, all the matter 
fields, transforming in the adjoint, get masses of the order of Mkk- 

The Yang-Mills 0-angle is holographically related to the integral of the RR one-form C\ 
over the circle S X4 . In [3] this field was taken as a probe of the original Z04-brane background. 
This approximation is fine if one is interested in keeping 9/N c very small and getting only 
the leading order corrections in this parameter. 

In this note we want to consider the full backreaction of the 0 angle on the background, 
with the aim to explore at deep how the IR Yang-Mills physics is affected by it. 

The relevant type IIA gravity action, in string frame, reads 


S = 


1 

2 *1 


d 10 xy/^g 


e ~W (n + 4 (defy 2 ) 




( 2 . 2 ) 


Here 2k q = (2n) 7 E s where l s = a fa! is the string length, F 4 = dCs is the RR four-form which 
is magnetically sourced by the N c Z04-branes, F 2 = dCj, and </> is the dilaton. 

The 0-backreacted background (see appendix [A] for a detailed re-derivation based 

on the Id reduced action) is a solution of the equations of motion derived from the action 
above. The string frame metric reads 



/«\ 3 / 2 

\r) 


\/ldo dx^dx^ + 



+ 



du 2 

T 


+ u 2 dQ 4 , 


where 


/ = !- + 

U A 


H 0 = 


1 - 


i4 0 2 
1 + 02 • 


(2.3) 


(2.4) 
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The background also includes a running dilaton, a RR one-form and a four-form field strength 
given by 

^ = 3 s ( 15 ) 3/4 Hq /4 , Ci = -^-dx A , F 4 = 3R 3 ui4 , (2.5) 

\K/ g s H 0 

with the flux quantization condition fixing the value of R as 



8r>if 9 ,AV R = {*g,N c yl%. 


( 2 . 6 ) 


Here 0 is a parameter which, as it will be clear in a moment, is proportional to 9. The 
original background at zero 9 angle, found in [T], can be readily obtained from the one above 
setting 0 = 0. Let us recall that in the expression above /x = 0,1, 2, 3 are the 1 + 3 Minkowski 
directions where the Yang-Mills theory is defined, d£l\ is the metric of a S A of radius one, u 
is the transverse radial coordinate u G [uo, 00 ) - holographically mapped to the RG scale of 
the dual field theory - X 4 is the compact coordinate of length 2nM^ K and R is a curvature 
radius. The isometry group of S' 4 is mapped into a global SO( 5) symmetry group in the 
dual field theory, which acts non-trivially on the Kaluza-Klein massive modes (signaling that 
these are, in fact, not decoupled in the limit we are considering). Finally g s is the string 
coupling and cu 4 is the volume form of the transverse S A , of volume V 54 = 87 t 2 /3. 

Notice also that, crucially, the S X4 circle shrinks to zero size when u = uq. Absence of 
conical singularities at u = u 0 is guaranteed if 


M 0 = 


4R 3 „ t2 1 
9 M '<K 1 + Q2’ 


(2.7) 


which thus implies that the (aq, u) subspace has the topology of a disk. 

Reliability of the background requires eC to be small: when this condition is violated we 
should better make use of the lid (“M-theory”) completion of the model. As it was shown 
in [TJ, in fact, the lOd solution can be obtained starting from an AdS-j x S 4 planar black 
hole solution of lid supergravity, reducing it on the M-theory circle and performing a double 
Wick rotation. We will make use of this picture in the following. 

As it will be reviewed below, the integral of C\ over S X4 at u —* 00 , is fixed to be 
proportional to 9 by the holographic dictionary. Then, since the S XA circle shrinks, a trivial 
solution like C\ rv_/ 9dx A is not allowed. Regularity requires in fact that C\ goes to zero at 
u = uo, a condition which is precisely satisfied by the solution above. 

The UV ’t Hooft coupling and the 9 angle of the gauge theory, can be related to the 
gravity parameters by considering the low energy limit of the D4-brane action 


S D a = -r 4 Tr 


d A x dx±e ^\J — det(G + T ) + r 4 


C*5 + -T 4 


C\ A Tr F A T , 


( 2 . 8 ) 


where T a Q = 2'Ka'F a p is proportional to the gauge field strength, t 4 = (27t) A l s 5 , C 5 is 
the electric five-form sourced by the branes (its field strength F 6 is the Hodge dual to F 4 ) 
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and G a p is the induced metric on the world-volume. Expanding the action to second order 
in derivatives, considering the UV asymptotics u —> oo, integrating over the compact X 4 

direction and comparing the resulting 4d action with the Minkowski version of (2.1), one 

461 


A — g N c — 2A4 — c 2gYf [ /[N c — ATrg s N c l s MxK 1 


9 + 2nk = — 


C u 


(2.9) 




where k is an integer. In the above expression we have introduced the parameters gvM and 
A 4 which are often referred to as the (UV) 4d gauge and’t Hooft coupling of the holographic 
model in the literature, despite the fact that they differ for a factor of 2 from the standard 
bare ones appearing in eq. (2.1). As we have observed above, A 4 ~ T S /M\ K parameterizes 
how far the holographic model is from pure Yang-Mills. 

The second relation in (2.9) defines 9 mod 2i r integer shifts (since the integral of C\ is 
gauge invariant only modulo 2 nZ) realizing the expected multi-branched feature of the dual 
gauge theory vacuum. Moreover, together with (2.5), it implies that the bare 9 angle is 
related to the 0 parameter of the background by 


0 = 


A4 

4:/r 2 


9 + 2kn 
N r 


A 

87T 2 


9 + 2kn 
N r 


( 2 . 10 ) 


As it is also suggested by eq. ( 2 . 1 ), in the’t Hooft limit, the corrections to the physics due to 
the 9 parameter, w.r.t. the 9 = 0 case, are actually weighed by the combination appearing 
in 0. Moreover, since this parameter depends on k, what we actually get on the gravity side 
is an infinite family of solutions corresponding to possible field theory vacua. 

As it was shown in [5] the curvature invariants of the background remain small if |@| <C 
A 4 / 4 . A similar constraint is obtained by studying the tension of the domain walls (identified 
with wrapped H 6 -branes) separating the various A:-vacua in the dual field theory: if |@| is too 
large, the (metastable) vacua become unstable. Our results will thus be reliable provided 
0 satisfies these bounds. Actually, as we will see in the following Section, on the field 
theory ground-state the possible values of | 0 | turn out to be bound after minimizing the 
multibranched energy density. 

Let us conclude this Section by recalling that there is an alternative way of getting the 
holographic relations (2.9). The one-instanton action 87 r 2 /g 2 + i9 in the gauge theory, is 
mapped into the Euclidean on-shell action for a /X)-bra lie wrapping the S X4 circle])] In 
Minkowski, the relevant DO-brane action is given by 


Sdo = -To / dx 4 e %+ To / Cl 


( 2 . 11 ) 


6 Remember that f TrF A F = (1/2) f d 4 xTrFF. 

7 This relation is readily obtained by considering the on-shell value of the Euclidean version of the D 4- 
brane action written above, on a 4d one-instanton solution F = F (the solution F = —F corresponds to 
an anti-DO-brane), with f d 4 xTiFF = 167t 2 . Indeed, the f C\ A Trf A F term gives one unit of Z?0-brane 
charge on-shell, since r 0 = (27r/ s ) 4 T4. 
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where r 0 = Ij 1 is the DO-brane charge. Using the u —>■ oo limit of the lOd background 
introduced above, going to the Euclidean frame, and performing the identification with the 
one-instanton action, one precisely gets the relations (2.9). 


3 Observables 


In this Section, using standard holographic methods, we will extract relevant information on 
the physics of the Yang-Mills theory dual to the 0-backreacted background given above. The 
main aim is to study the 9 dependence of interesting physical observables, with an eye to the 
available results for pure Yang-Mills on the lattice (see e.g. 0 ). As we have already stated 
in the Introduction, studying the physics at finite 6 is challenging on Euclidean lattices due 
to the fact that, as shown by eq. (2.1), the 6 term in the Euclidean Lagrangian is imaginary. 
Lattice results are thus obtained either computing the coefficients in series expansions around 
6 = 0 or by extrapolations from imaginary 9 angle, again around 9 = 0. Correspondingly 
only the first few terms in powers of 9 are generically computable]^ 


In the holographic model, instead, in the limits where the solution is reliable, the physics 
can be easily studied and the results exactly given at any order in 0. Focusing on the small 9 
regime where a comparison with lattice results is sensible, we will see that the physics of the 
holographic Yang-Mills model precisely matches, at least qualitatively, with that of lattice 
Yang-Mills. This coincidence let us try to adopt the holographic model as a way to predict 
the behavior of the subleading order terms in the 9 expansion, for the realistic model. 

Moreover, for certain quantities, like the entanglement entropy, which cannot be computed 
yet on the lattice, not even at 9 = 0, the results we obtain for the holographic model are 
unique and could hopefully provide useful benchmarks for pure Yang-Mills. 


3.1 The ground-state energy 

One of the basic entries in the holographic dictionary is the relation between the field theory 
partition function and the (renormalized) on-shell gravity action. The ground-state energy 
density f{9 ) of the Yang-Mills theory can thus be obtained through the relation (valid at 
large N c and A 4 ) 

Z = e~ nf{6) « e" s sTn- sh eii ; ( 3 . 1 ) 

where V 4 is the (infinite) 4d Euclidean spacetime volume and the renormalized on-shell 
Euclidean gravity action is given by (see also [18] for a recent account of the holographic 
renormalization in the model at hand) 

“-h?on—shell = $E + S GH + S c _ t . . (3.2) 


Effects of the 0 term in certain Yang-Mills theories can be studied with semi-classical methods m- 





Here 


S R = 


2 kl L . 


d w xyjg 


e~ 2 * (R + 4(d<P) 2 ) -~\F,\ 


-1\f 2 I 


is the Euclidean version of the action (2.2|). Moreover 


S GH = -^ I d 9 xVhe~ 2<i> K , 


(3.3) 


(3.4) 


is the Gibbons-Hawking term, where h is the determinant of the boundary metric (the slice 
of the lOd metric (2.3) at fixed u = u\ with u\ —> oo being the radial position of the UV 
boundary) and K is the trace of the extrinsic curvature of the boundary 


K = 

V9 


V9 \ 

yjQuu J 


u=u\ • 


(3.5) 


Finally, the counter-term action is given by [19 


qbulk _ x 

D c.t. — T/2 


1/3 ' 

G_ 

Ji 


d 9 xVh-e~ 7 ^\ 


(3.6) 


Evaluating the above terms on the 0-backreacted background given above, one discovers that 
there are no ^-dependent divergent terms and that the final result reads formally as that for 


0 = 0 


s ren 


E on—shell 


2fc o g\ 


:V 4 


M kk Vs4 2 


(3.7) 


The 0-dependence comes out explicitly through the relation (2.7). Using the holographic 
map (3.1) and expressing everything in terms of field theory quantities, from (3.7) it follows 
that the ground-state energy density of the Yang-Mills theory is formally given by (see also 
0 ) 

/(e) = 0.8) 


3¥ (i + e 2 ) 3 ’ 


where the negative sign is in agreement with the negative Casimir energy one expects to 
arise from the compactihcation of a supersymmetric theory with supersymmetry breaking 
boundary conditions [20]. 

Actually, since 0 is proportional to 0 + 2kn, for a given value of 0 the true vacuum energy 
is obtained by minimizing the previous expression over k 


/(0) = min fc /(0). (3.9) 

As a result, the ground state energy density (see figure [Tj) turns out to be a periodic function 
of 0, as expected. To any given interval, of length 27T, of possible values of 0, it will correspond 
a precise value of k. For example, k — 0 when 0 e (—7r, 7 t), k — 1 when 0 6 (— 3n, — n), 


9 



















f(9) ~ 1(0) 

l/(0)j 




-471 


-3tt 


-271 


2n 


3n 


4n 


Figure 1: The (normalized) ground state energy density at finite 8. Different colors (online) 
correspond to different branches of the vacuum energy: k = 0 when 8 E (— ir, 7r), k = 1 when 
8 E (—37r, —7r), k = — 1 when 8 E (7r, 37 t) and so on. At 8 = ±7r and odd multiples of these values, 
a CP breaking first order phase transition separates the different fc-branches. The plot has been 
obtained from the function in the text setting A4/(47r 2 lV c ) = 0.1. 


k — — 1 when 9 E (7r, 37 t) and so on. All in all, in the ground-state 1 6 + 2kn\ < it, so that 
|©| < A 4 /(47riV c ). 

The theory experiences a first order CP-breaking phase transition when passing from one 
branch to another. The transitions happen at 6 = ±7r and odd multiples of these values: CP 
symmetry, which would be preserved at these points, breaks spontaneously by the choice of 
a branch. In the small 9 limit (which corresponds to the k — 0 branch) one gets 


m - m = nXgO' 2 


9 


9 4 


1 + + 0(9 


with the topological susceptibility given by [Ij 

Xg = 

and the expansion coefficients given by 

X 2 

h = ± 

2 8tt 4 ’ 


‘ N? 


4(37t) 6 ’ 


6,= 5A 1 


n; 


384ir 8 


(3.10) 


(3.11) 


(3.12) 
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The & 2 n coefficients provide relevant physical information as they are related to the zero- 
momentum n-point connected correlation functions of the topological charge density at 9 = 0 
[2], These are also phenomenologically interesting quantities: for instance 62 (for N c = 
3) is related to the 7 / — 7 / elastic scattering amplitude. In turn, the b 2n coefficients give 
the moments of the probability distribution of the topological charge Q. Non-vanishing 
coefficients imply that this distribution departs from a simple Gaussian one. 

It is interesting to compare these results with those obtained for pure Yang-Mills on the 
lattice (see [[2, [21]). Calling b 2n (unbarred) the coefficients multiplying just the 6 2n powers, 
the following values have been found for N c = 3,4, 6 

62 (lattice) = -0.026(3), -0.013(7), -0.008(4), (3.13) 

yielding a large N c estimate for the barred coefficient 

b 2 «-0.2. (3.14) 

It is remarkable to notice that the sign of this coefficient precisely matches with that obtained 
in the holographic model. 

Concerning the value of 64 on the lattice, at present the errors are such that its sign 
cannot be safely determined: recent results [ 8 ] just provide a bound on its absolute value 
| 7 > 4 1 < 0.001. Provided the qualitative matching of the holographic model and the pure Yang- 
Mills theory persists at subleading order, our results suggest that 64 should have a positive 
sign. 

3.2 Rectangular Wilson loop: the string tension 

As it was pointed out in [22], in the large N c , large’t Hooft coupling regime the VEV for a 
Wilson loop on a contour C is holographically given as 

(W[C}) ~ e~ s ™ ; (3.15) 

where S^g is the renormalized Nambu-Goto on-shell action for a fundamental open string 
whose end-points span the contour C. For a rectangular contour with sides of length T along 
the x° = t direction and length L along one space direction, say x = x 1 , computing the 
Wilson loop allows to obtain the static quark-antiquark potential V(L). In the large L limit 
in fact 

{W[C}) « e“ Ty(L) . (3.16) 

Linear confinement of the chromoelectric flux tube implies that V(L ) = T S L at large L : the 
Wilson loop has an area law and T s is the string tension. 

The holographic computation of the rectangular Wilson loop proceeds as follows. Since 
we are interested in comparing the result with that obtained in pure Yang-Mills on the 
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lattice, we just consider a Wilson loop carrying no spurious global symmetry charges: the 
corresponding open string will thus be point-like in the transverse S' 4 as well as in S XA . 
The relevant string embedding can thus be chosen, in the static gauge, as r = x° G [0, T], 
cr = x e [—L/2, L/2], u = u(x) so that the Nambu-Goto action is given by 


Sng = 


2i:a' 


clTcllJy/ fjrrfjrra 


2 not' 


-T 



-goo(9xx + 9uuU'{x) 2 ) , (3.17) 


where g is the (induced) string frame metric. The Euler-Lagrange equations for u(x) give 
the actual string prohle to be used for computing the on-shell action. As it has been largely 
discussed in the literature, for the kind of background we are considering, in the large L limit 
the string prohle is bath-tube-shaped: in order to minimize its energy the string wants to 
stretch as much as possible over the u = Uq region where its effective tension, proportional 
to V~9oo9xx, is minimal. Then, around the extrema x = —L/2 and x = L/2, the string 
will move approximately vertically up to the UV cutoff u = u\ —)• oo where it is attached to 
a probe brane. The two vertical branches of the string give divergent contributions to the 
action when the cutoff is sent to infinity. These contributions can be interpreted as due to 
the (infinite) masses of the static quark-antiquark pair and must be subtracted in order to 
get the renormalized Nambu-Goto action. All in all, in the large L limit one gets 


5' 


NG 


27red 


goo9xx | u=u o L , 


(3.18) 


which, using the holographic map (3.15), implies that the Wilson loop obeys an area law, 
with a string tension given by 


T„ = 


‘I'KCt! 


V 9oo9xx |u=uo • 


On the ^-deformed background (2.3) this reai 

_ 2A 4 2 

s 27ti Mkk 


(1 + 0 2 ) 2 


(3.19) 


(3.21) 


As it has been shown in the previous Section, one has to read this expression carefully: 
actually O = 6 + 2kn and, in the ground-state, for a given value of 9 in an interval of length 
27r one has to fix the corresponding value of k. 


Setting k 


0 and considering the 6 —> 0 limit one gets 


T, = 


2\ 


27tt 


4 m 2 kk 


Ai 


e 2 


8t r 4 A" 2 


+ 


3A 4 


e 4 


256tt 8 N c 4 


+ 0(9 6 


(3.22) 


9 In [53] the leading stringy corrections to this result have been derived. Following the rule for the 0 
corrections mentioned in the Introduction, the full results would then read 


T„ = 


2A 4 

27V 


M 


KK 


(1 + 0 2 ) 2 




0 2 ) log 2 


(3.20) 
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It is remarkable that the 0(9' 2 ) correction is negative. This, together with the expected 
scaling with 9/N c , is precisely what has been found on the lattice in [2] 


T s iat — T s i at ( 0 ) 


9 2 

1 + S 2^ + --- 


where, using the large N c extrapolation of the N c = 3,..., 6 results one geti 

s 2 ~ —0.9. 


10 


(3.23) 


(3.24) 


Moreover, using the holographic result (3.22) we are led to predict that the 0(9 4 ) correction 
to the string tension in pure Yang-Mills should be positive. 


3.3 ’t Hooft loop and oblique confinement 


The rectangular ’t Hooft loop computes the monopole-anti-monopole potential. At 9 = 0 
one expects that confinement of the chromoelectric flux is associated with screening of the 
magnetic charges. The situation is expected to change, instead, when the 9 angle is turned 
on. Let us see what the holographic results tell us about this. 

The “chromomagnetic string” in the model we are considering is a probe T>2-brane wrapped 
on S X4 [241 [IT] . In order to follow standard conventions, it is useful to write the relevant 
part of its action as 

S D 2 = ~T 2 [ d 3 £e- < V-det(3 + J r ) + T 2 [ C { A J= , (3.25) 


where T 2 = g~ l r 2 = (2 / n)~ 2 g~ 1 l~ 3 is the L>2-brane tension, e~^ = g s e~^, Ci = g s C\ and 
T = 2na'F is a U( 1) gauge field on the brane world-volume. Notice that due to the non¬ 
trivial Ci potential in the background, the equations of motion for F are not solved by the 
trivial solution F — 0. Instead of working out these equations, it is easier to make use of the 
relation between the action above and that of a wrapped M2-brane in the lid completion 
of the background (as we will see, only the low energy details of the latter play a role in the 
computations, hence working with the M-theory completion is equivalent to working with 
the lOd one). This relation can be obtained using the following procedure (see e.g. [25]). 

One introduces an auxiliary vector field t a on the H2-brane world-volume and rewrites 
the action as 


s D 2 = -t 2 / d 3 e 


det(g a/ 3 + e 2 H a tp) - -e Q/37 (C' Q - t a )F^ 


(3.26) 


Integrating out t a one gets back the original action. If instead we treat the gauge field A in 
T = dA as a Lagrange multiplier we see, from its equation of motion, that 


e^dpiC^-ty) = 0 , 


(3.27) 


10 Calling S 2 the coefficient multiplying just 6 2 , on the lattice on finds, for N c = 3, S 2 = —0.08(1). 
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which implies that C a — t a = d a y with y being a scalar. Using this equation one discovers 
that the 142-brane action becomes equivalent to the M2-brane action 

S M 2 = - ( 27r )2/3 J det G , (3.28) 

1 /3 

where In = l s g s and G is the pullback on the brane world-volume of the lid metric 

ds\ x = e~^ds \o + ea \dy — C XA dxfi) 2 . (3.29) 


Here dsf 0 is the lOd string frame metric given by (2.3) and y is identified with the compact 
eleventh direction of M-theory of length 2nR y = 2n g s l s , compactifying over which one gets 
back the lOd IIA string theory background. In lid the whole background reduces to the 
metric above and a four form held strength, which solve the equations of motion of the lid 
supergravity action. 


As already pointed out in H, there are two relevant cycles one can identify in the (aq, y) 
subspace. The first one is defined by the equation x 4 = const: this is the cycle over which 
one reduces M-theory to get the type IIA model. An M2-brane wrapped over this cycle 
reduces to a fundamental string. The second cycle is what we are interested in: it is the 
minimal volume contractible cycle over which we have to wrap the M2 -brane to get the 
chromomagnetic string. Looking at the UV (u —> oo) asymptotics of the background, we see 
that this cycle is defined by 

y = @x 4 + const . (3.30) 


The metric on this cycle reads, thus 


ds 2 


.2 1 4 

3 y44 + e 3 


0 0 2 



dx 4 . 


(3.31) 


Notice that the volume of this cycle is minimized at u = uo, where it takes a finite value if 0 
is different from zero. In this case, thus, the wrapped M2-brane action, in the limit where 
the monopole-antimonopole distance L is very large, is minimized by a bath-tube shaped 
configuration, precisely as in the case of the Wilson loop fundamental string. Using the 
same logic as in that case, one gets, for the chromomagnetic string tension, the following 
remarkable expression 


T 

1 rr 


1 n 4~2 x \d + 2k7T\ _ ^ \9 + 2kn\ 
27 ^ MkkA \i + Q^- 1s 2n 


(3.32) 


Hence, at any order in 0 within our approximations, there is a very simple relation between 
the chromomagnetic tension T m and the string tension T s . Notice that at 9 = 0 {k = 0) this 
relation implies that T m = 0: as expected there is no area law for the ’t Hooft loop and the 
magnetic monopoles are screened. At finite 9 , instead, the ’t Hooft loop generically shows 
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an area law. The objects which are screened are actually dyons, particles of electric charge 


—p and magnetic charge q. Using the above relation we see that the “string” tension for the 
dyons is (on the k — 0 branch and taking 6 > 0 ) 



(3.33) 


so that dyons are screened provided that 6 = 2Tr(p/q)\ this is precisely the so-called oblique 
confinement. 

3.4 The scalar glueball mass 

Confinement in Witten’s holographic model is accompanied, just as in pure Yang-Mills, by 
the occurrence of a mass gap in the (glueball) spectrum. In pure Yang-Mills the lightest 
glueball is a C and P even scalar mode 0 ++ . The 6 dependence of its mass has been studied 
on the lattice in [ 6 ]: there, the 0 ++ mass was found to decrease quadratically with 6 around 
9 = 0. As customary, the regime investigated by holographic means is different from the 
lattice one. Nevertheless, we are after the qualitative behavior of the observables. This is 
hoped to be a stable feature of the theory. So, let us compare the lattice result with the 
holographic model. 

At 6 = 0 the spectrum of glueballs in Witten’s model has been studied in many papers 
(for example in [20, [27, 20, I29j). There are actually two light 0 ++ modes. The lightest 
one was dubbed “exotic” in m, as it comes from a metric perturbation involving, among 
many others, the compact £4 direction of the background; it is sometimes called “mode S” 
in the literature. The second light mode was dubbed “dilatonic” in [3D], as it involves a 
fluctuation of the dilaton in the ten dimensional geometry; it is customarily called “mode 
T”. It is degenerate in mass with the tensorial 2 ++ and vectorial 1 ++ glueballs, due to the 
symmetries of the background - it comes from the extra flat Minkowski direction in the 
eleven dimensional origin of the background, so the equation governing its mass is precisely 
the same as the one for the Minkowskian tensorial perturbations. 

Since it is lighter, not degenerate with the 2 ++ and common to some deformations of 
Witten’s model [29], at first sight it would seem that the “exotic” mode S is to be considered 
in the comparison with Yang-Mills data. This could actually not be the case. As a first point, 
from inspection of the DBI action of a DA- brane, one can see that the operator TrF 2 sources 
both the S and T 0 ++ modes. Furthermore, the “exotic” polarization could not survive in 
the limit where one decouples the KK modes, as it comes primarily from the excitations 
along the KK direction £ 4 . Moreover, the extrapolation of lattice data at large N c is still 
imprecise [29]. Finally, a recent study of glueball decays has shown that the “exotic” mode is 
too broad (and too light) to be compatible with the scalar glueball candidates in QCD [3D]p] 

11 See also 31 1 . 
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All in all, although the situation is not firmly settled yet, it seems possible, if not likely, that 
it is the “dilatonic” T mode and not the “exotic” S mode the correct one to be compared 
with QCD - it should be the only one sourced by TrF 2 in the KK mode decoupling limit. 
For these reasons, we will concentrate on the “dilatonic” T mode in the following. 


It is convenient to work in eleven dimensions, so that the only field involved in the dis¬ 
cussion is the metric. The 0-backreacted background it is given by eq. (3.29). The general 
equation for the metric fluctuations (also at finite 6), obtained by linearizing the lid Einstein 
equations on the background, reads 


1 

2 


V„V^ + -V 2 /r ab 


V c V (a V 


6 

4 


9-ab 


0. 


(3.34) 


Here a, b, c, ... are seven dimensional indexes referring to the reduction of the eleven di¬ 


mensional background (3.29) on the four sphere, while //, zz, ... = 0, ...,3 refer to the four 


Minkowski directions. The fluctuations can be factorized as h a b = H a b{u)e lkx , with k 2 = 
—Ad 2 ] we work in the frame where k^ = cadf. 

The “dilatonic” mode T at 6 = 0 corresponds to a traceless excitation H a b = e a bH(u ) with 
diagonal components only along the spatial directions with /i = 1, 2, 3 and r/j^] The key 
point in our discussion is that it is sufficient to keep the traceless condition also at finite 6 


in order for the mode to be compatible with equation (3.34), provided one single equation 


for H{u ), giving the mass spectrum of the glueballs, is satisfied. The correct ansatz for the 
fluctuation is 

Hab(u) = ^-H(u) diag ( 0,1,1,1,0,- 3 2 ,0 ) , (3.35) 


R 


i + e 2: 


giving the equation 


H"{u) + 


4 u 3 


Un 


-H'\ 


M 2 R 3 


u 


u(u 


u n 




Un 


H(u ) = 0. 


(3.36) 


This has acceptable solutions (regular at u = uq and normalizable at u —>• oo) only if M 2 > 0 
and the spectrum is discrete. This is the holographic realization of the mass gap pQ. 


The fact that the ansatz (3.35) is sufficient was not a-priori guaranteed, since the metric 


at finite 6 has non-diagonal entries which are expected to cause a mixing among the 6 = 0 
modes. Indeed, the fact that generically the mixing happens was found in [32j. There, the 
backreaction of the RR one-form on the metric was not considered, but taking into account 
its 6 dependent background value it was shown that the gravity mode dual to the 0 _+ glueball 
mixes with the T and S modes (while the latter obviously satisfy the 9 = 0 equation, as the 
metric was unchanged). The mixing between the 0 ++ and 0 h glueballs has to be expected, 
as 9 breaks CP invariance. In fact, somewhat surprisingly, in [32] it was found that, despite 


12 At 0 = 0 the equation of motion for this fluctuation coincides with that of a minimally coupled scalar. 
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the mixing of the modes, the masses of these glueballs are not changed if the backreaction 
of F 2 on the geometry is not taken into account. 

What we have found here is that, despite the full backreaction of 9 on the geometry, the 


small modification of the ansatz for the mode T in (3.35) still preserves its equation: in fact 


(3.36) is identical to the equation at 9 = 0. This fact was observed in [5] for the lightest 2 ++ 
mode (although in that case no mixing is actually expected). Note that the 2 ++ and 0 ++ T 
modes are still degenerate in mass at finite 9. 

The fact that the equation for a mode is unchanged by 9 does not imply that the mass 
of the glueball does not depend on 9. In fact, the masses are given in units of Uq/R 2 , 


whose relation with the physical scale Mkk is 9 dependent: from equation (2.7) it follows 


immediately that the 9 correction to the mass of the lightest relevant glueball is given by 

M(0 = 0) 


M(0) = 


(3.37) 


VT+W ' 

As usual this expression implies that M(9) is a periodic function. Focusing on the k — 0 
branch and expanding around 9 = 0 one thus finds 


M(9) = M{9 = 0 ) 1 - 


Al 


9 2 3A 4 9 4 


+ 0(9 6 


(3.38) 


327T 4 N 2 ' 20487T 8 IV 4 

The leading correction in 9 2 has a negative sign, precisely as it was found in lattice Yang-Mills 
[ 6 ], where, for N c = 3 

M lat (9) = M lat ( 0) (1 + g 2 9 2 + 0(9 4 )) , g 2 = -0.06(2). (3.39) 

The holographic model predicts that the 0(9 4 ) correction has a positive sign. Again, it 
would be nice to check this expectation on the lattice. 


3.5 The mass of the baryon vertex 

In Witten’s holographic model, a baryon vertex is identified with a iA4-brane wrapped on 
S ' 4 and localized at the radial position corresponding to the deep IR of the dual held theory 
[33]. Using the 0-backreacted background, we can easily study how a finite 9 term affects 
the mass of the baryon vertex. The wrapped ZM-brane action reads 


Sjj 4 = —t 4 / dx°dtt 4 e ^\J det g§\ u =u 0 = —"is / dx° 


(3.40) 


where 74 = (27r) 4 l s 0 is the i74-brane charge. Using formulas (2.3)-(2.6) in (3.40), integrating 


on S ' 4 and substituting the expressions in (2.7), (2.9), one finds that the baryon vertex mass 
is given by 


m B (Q) = —^ N c M K k 


m B ( 0) 


277r - + 02 ) 3/2 (! + 02 ) 3/2 > 

which is a periodic function of 9. Thus, the mass of the baryon vertex decreases with 9 


(3.41) 
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3.6 Finite temperature: the confinement-deconfinement transi¬ 
tion 


Going to finite temperature in Witten’s holographic Yang-Mills model, amounts to compact- 
ifying the Euclidean time on a circle of length /3 — 1/T. As in the 9 = 0 case, there are two 
possible solutions of the gravity equations of motion for which this condition can be satisfied 
(see e.g. [34]). 

One solution has precisely the same structure of the one discussed above, the only difference 
being just in the fact that its Euclidean continuation has a compact time circle. This solution 
corresponds to the confined phase of the theory. The free energy density (hence minus the 
pressure) in this phase has the same expression as in the T = 0 case (3.8) 


2 AYA 4 M$ <k _ /(0) 

3 7 vr 2 (1 + 0 2 ) 3 (1 + 0 2 ) 3 ' 


(3.42) 


The usual gravity solution corresponding to the deconhned phase has a black hole event 
horizon and the x 4 circle does not shrink anymore J 73 ] This allows F 2 = 0, i.e. C\ = 9dx 4 
to be a solution of the background equations of motion. As a result the background metric 
does not depend on 9 and it is given by 


/ u \ 3 / 2 


/ u \ _3 / 2 

du 2 0 0 

y 

f(u)dxl + dx a dx Ll + dx'l 

+y 

—-b u 2 d{}j 

L/(«) J 


~ 7 /^ 

/(«) = 1 - -i ■ (3-43) 

u 6 

The dilaton and the F 4 form do not change w.r.t. the T = 0 solution. In the equation 
above, a = 1,2, 3. The (x 0 , u ) subspace is topologically a disk - with the S XQ circle smoothly 
shrinking to zero size at u = Ut (the position of the horizon) - provided we identify 


9 (3 2 ut = I67r 2 i? 3 . 


(3.44) 


The fact that this background actually corresponds to a deconhned phase of the dual held 
theory is readily seen recalling, for example, what we got for the Wilson loop: the would-be 
string tension in this case would be zero since goo9xx\u=u T = 0. 

Using the Bekenstein-Hawking formula one gets the black hole entropy density (which is 
holographically mapped to that of the dual held theory) 


256Y c 2 7t 4 A 4 5 
729 M] (K 


(3.45) 


13 In j25] it has been proposed an alternative solution for the deconhned phase and [3B] presents a study 
of instantons on that background. 
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from which we can deduce the free energy density (hence minus the pressure) using s = 

-( df/dT) 


fdec Pdec „ 


1 256fV?7r 4 A4 


6 729 M 2 kk 


rpt 


(3.46) 


Notice the peculiar scaling with T, which follows from the higher dimensional UV completion 
of the model. Since the free energy is independent from 9 , one gets that in the holographic 
model at hand, the topological susceptibility is zero in the deconhned phase. This is actually 
what one expects in the large N c limit of Yang-Mills, where the susceptibility is exponentially 
suppressed] 44 ] 

Comparing the free energies (or pressures) in the two allowed phases on gets that the dual 
field theory experiences a first order confinement-deconfinement phase transition at a critical 
temperature T c obtained by solving the equation p = pdec{T c ) where the pressure p in the 
confined phase is the one given in (3.42). The result is 


T C (Q) 


MkK 1 
27T y/l+W ' 


(3.47) 


As usual we have to recall that 0 is actually proportional to 9 + 2kn and that, as we have 
observed studying the vacuum energy density, a given interval (of length 27 t) of possible values 
of 9 is related to a precise value of k. As a result, plotting the actual critical temperature 
T c (9) one gets a periodic behavior like that shown in figure [2] It is remarkable that the 
holographic model explicitly realizes the expected but still conjectural (see e.g. [7j) “arcade” 
structure, allowing, in turn, to study the behavior of the phase diagram near the critical 
points (which turn out to be triple points) where the first order deconfinement transition 
line meets the first order lines related to the CP breaking transitions in the confined phase P’1 
It is tempting to conjecture that the same cusped structure should appear in pure Yang-Mills 
at large N c . 

Focusing on the k = 0 branch and expanding around 9 = 0 we find (at 0(9 2 ) the result 
was already found in [18]) 


U9) 


Mkk 

2tt 


At 9 2 3A 4 9 4 

32^0 AY + 2048vr 8 Nj 


+ 0(9 6 ) , 


(3.48) 


14 In the deconfined high temperature phase the topological sector of the Yang-Mills theory is well described 
in terms of an instanton gas. The same description does not hold anymore in the confined phase. This is 
encoded in the dual holographic picture (see e.g. mo 13E] and references therein). As we have seen, instantons 
are identified with DO-branes wrapping the S X4 cycle. In the deconfined phase, this cycle does not shrink 
and the -DO-brane configuration is stable. In the confined phase, instead, the S Xi circle shrinks and the same 
configuration is unstable. 

15 Let us focus for example on that at 6 = 7r. It is easy to show that the phase diagram at that point 
has a cusp: when 8 —> 7r =l= , T c (9)/T c ( 0) ss (T c (7r)/T c (0)) + A(tt ± 8 ), where A > 0 is given in terms of the 
parameters of the model. 
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T c (9)/T c (0) 



Figure 2: The (T, 6 ) phase diagram. Different colors (online) correspond to different branches of 
the vacuum energy: k = 0 when 6 £ (— 7r, 7 r), k = 1 when 0 £ (— 37T, — 7 t), fc = —1 when 0 £ (7r, 3n) 
and so on. The critical temperature has cusps at 6 = ± 7 r and odd multiples of these values. The 
cusps are actually triple points where the deconfinement first order transition line meets the CP 
breaking first order transition dashed lines separating the different /c-branches in the confined phase. 


where, again, both the quadratic dependence and the sign of the leading 0 -dependent cor¬ 
rection, agree with what has been found on the lattice [7j. There, for N c = 3 it has been 
found that 

Tc(e) lat = T c (0)m [1 - ReO 2 + 0(0 4 )] , R, = 0 . 0175 ( 7 ). ( 3 . 49 ) 

Our model, then, predicts that the 0(9 4 ) correction has a positive coefficient. It would be 
nice to check this prediction on the lattice. 

3.7 Entanglement entropy 

In this section we compute the dependence of the entanglement entropy on the 9 angle, as 
a tool to characterize the phases of the theory. As customary, we take the entanglement 
entropy as a measure of entanglement between two physically disjoint spatial (t = const.) 
regions in the boundary theory: region A and its complement, region B. These regions are 
separated by a given surface 7 , whose shape can enter the result for the entanglement entropy. 
We will mostly consider the “slab” and the ball geometries for region A. The results will be 
that in the UV the entanglement entropy is unaffected by the 9 angle, while it is suppressed 
in the IR; the critical length at which this confining model exhibits a first order transition 
in the entanglement entropy increases with 9. 

The holographic prescription for the entanglement entropy [32] (see also mm consists 
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in calculating the (Einstein frame) area of the minimal (eight dimensional in our setup) bulk 
surface T which is anchored to the given surface 7 at the boundary separating regions A and 
B. This is achieved by minimizing the action (which we write in terms of the string frame 
metric (Jmn = e ^ 2 9MN, since this is the one we have worked with in the previous Sections) 

S = —[ d 8 cre _29 V- det(Gi n d), (3.50) 

tv J 


where G^ — 87 r 6 g;?Z® is the ten dimensional Newton constant and G in d is the induced string 
frame metric on the eight dimensional (fixed time) bulk surface T. 

In our background it is straightforward to realize that the form of the metric and dilaton 
(2.3), (2.5) are such that the factors of H 0 , carrying the explicit dependence on 0, cancel 
out from the expression \J — det(G; n d) for any surface T. 

This can be seen also by performing a dimensional reduction (over the compact transverse 
spaces S ' 4 and S X4 ) to five dimensions, where the metric does not depend on H 0 anymore. 
The canonically normalized Einstein frame 5d metric reads, in fact 


ds%E = /(«) 1/3 Q) 


U\5/3 


(IXi/llf + — 


R\ 6 du 2 


u 


/(«) 


(3.51) 


with the 5d Newton constant defined as = G N ] R 4 V 54 2ttM k x k . The holographic entan¬ 

glement entropy can be thus equivalently deduced by minimizing the area 


S(5) = 


4 G 


N(5) 


d ^ \J det by 5 jind , 


(3.52) 


where G-'(, 5 )i n d is the pullback of the 5d Einstein frame metric (3.51) on the surface. 

All in all, the dependence on 0 just follows from the explicit expression of the stringy 


variables in terms of held theory ones, formula (2.7). 


As explicit examples we consider the “slab” geometry, for which the entanglement entropy 
has been derived in HU, and the ball geometry considered in ra- We can just repeat the 
analysis in those papers, to which we refer the reader for details, with minimal modihca- 
tionsE3 


3.7.1 The “slab” geometry 

This is the simplest geometry to consider: region A is the direct product of M 2 times an 
interval of length l EH]. The minimal surface in the bulk with these boundary conditions 
has distinct features at small and large l. For small l it is determined by a connected surface 
extending in the bulk up to a minimal radial position u * > Uq. For large l the minimal 

16 Note that the A = g s N c in HUE] is related to our A 4 by an inverse 2ttMkk factor. 
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surface is given by two disconnected pieces, anchored at the boundary to the two extrema of 
the segment and extending in the bulk all the way up to Mo- The critical length l c at which 
the (first order) transition between these two configurations happens could be viewed as a 
probe of the scale of confinement, since the connected surface explores a UV (“deconhned”) 
region, while the connected one a IR (“confined”) region. 

The entanglement entropy corresponding to the disconnected solution, dominating at large 
l, can be derived analytically and reads 


Sdis = 


Vo 


Vo 


4Gv(5)-R 


. u. 


Un = 


4 G 


-R 3 V s 4 


2tt 
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Mi 


l M, 


Un 


KI< 


(3.53) 


where V 2 is the (infinite) two dimensional volume of M 2 and u Q0 is a UV cut-off. In terms of 
held theory variables it can be written in the following form 


S dis = v 2 


A 4 -V 2 


7 t 3 5 (1 + @ 2 ) 2 


Mkk 



(3.54) 


where e = Uq/uoc —>■ 0 accounts for the UV divergence. Sdis has the typical parametric 
dependence of a 4d theory, which is the case in the IR of the Witten model. Also, the result 
does not depend on /, so that dSdi S /dl ~ 0(N®), while it is very sensitive to the 6 angle, 
which affects the IR of the theory. Notice that, as usual, the previous expression has to 
be read by taking into account the 6, k dependence of 0. As a result, Sdi s (0 ) is a periodic 
function. 


The connected surface allows to get an analytic form only for small /; the UV divergent 
part of the entanglement entropy is precisely as in (3.53), while its finite part reads 


5, 


conn, fin 


2 6 7 T 3 / 2 r [|] 5 A 4 iV 2 

2 3 r[A]5 (/ M KK y 


Ml K - 


(3.55) 


The dependence on l is characteristic of a 6d theory, consistently with the UV of Witten’s 
theoryp] Note that in this case dS conn /dl rv_/ o{nd [nj. Moreover, there is no dependence 
on 6, consistently with the fact that the UV of the theory is not sensitive to the shift in the 
vacuum energy. 

This fact provides immediately the behavior of the critical length l c with 0. In fact, 
comparing the two types of solutions, eqs. (3.54) and (3.55), it follows immediately that 

l c = Z c ,oVl + 0 2 , (3.56) 


being l C) o ~ 1.288/Mkk the critical length at 0 = 0. This expression is confirmed by 
the numeric comparison of the entanglement entropy of the disconnected solution with the 
connected one (l is not forced to be small in this case), see also figure [3] (showing the swallow- 
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Figure 3: The difference between the entanglement entropies corresponding to the connected and 
disconnected solutions as a function of the length l, for 0 = 0,0.7,1.4, 2.1 (left to right). Notice 
the typical swallow-tail behavior. 


tail trend typical of first order transition^]) . So, the critical scale for the phase transition, 
which is a probe of confinement, behaves as I c Mkk ~ \/l + @ 2 . The usual considerations on 
the k, 9 correlated values, imply that the actual critical length l c (Q) is a periodic function. 

The main effect of the 9 angle is to enhance the vacuum energy (in the “confined” phase). 
As such, the length at which the UV (“deconfined”) regime is probed is (in each branch) 
larger at larger \9\, as signaled by formula (3.56). 


3.7.2 The ball geometry 

In (2 + l)d the entanglement entropy of a large disc with its complement has a part inde¬ 
pendent of the disc radius r, which provides the so-called “topological entropy”, a measure 
of the topological order of a gapped theory. The natural generalization to (3 + l)d theories 
of this quantity, which in this case has no topological nature, is encoded in the entangle¬ 
ment entropy of a ball (with a two-sphere boundary) with its complement. This has been 
investigated for the Witten theory in [T2|. We can follow the latter reference to obtain the 
dependence of the entanglement entropy on the 9 angle. 

To begin with, it is interesting to note that the structure of the divergent terms in the 
entanglement entropy is richer than the one of the slab geometry. The divergent piece 

17 In pT] it is pointed out that using the fact that g s l s is the size R y of the M-theory circle, one can rewrite 
A 4 N^/Mkk ~ gNing the correct behavior of a 6 d theory compactified on a circle of length 2nR y . 

18 A similar behavior, for another non-local observable as the rectangular Wilson loop, has been found in 
QCD-like holographic models with massive flavors, see e.g. [421 . 
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exhibits both power-like and logarithmic terms 

(2vr) 4 R 3 


Sdiv 


1 2 2 3 3 1 -R 6 Uq 

U oo - “oo + log- 


8 r 2 


ZLr> 


3G^M KK 

In terms of field theory variables this can be rewritten as 

" 2 27 (1 + @ 2 ) 81 (1 + @ 2 ) 2 


(3.57) 


Srlii: 


Mrk , a 4 y 2 


35 (1 + 02)2- 


4e M\ k 


+ 


64 r 2 M 4 A - 


loge 


(3.58) 


where e = Uo/Wo —> 0 accounts for the short distance divergence. 

As in the previous section, there are two competing solutions with different topology and 
a phase transition between the two at a certain critical radius. The large r behavior is 
dominated by a cylinder-like solution, for which the entanglement entropy reads 


20 


S C yl — 


(2tt) 4 R 3 

3GnMrk 


1 2/ 2 2\ 1 0 

2 r Kc - «o) + ~Y~ 


n'U'OC 0 7T 

-6-b 2 + — + log 3 

u 0 V3 


(3.59) 


The first, r 2 -dependent part is the analogous to the result for the disconnected surface in 
the slab geometry (3.54): the two expressions coincid^j if we denote the volume of the two- 
sphere as V 2 . But in the present geometry there are also subleading terms. In particular, 
the r 2 -independent term in (3.59), which in (2 + l)d is by definition the (opposite of the) 
topological entropy, reads 


S r 2 —imrl 


A 4 1V 2 


54 (1 + © 2 ) 


6 7r 

__ + 2 + _ + log3 ^ 


(3.60) 


Again, since this is a IR quantity, it is very sensitive to the presence of the 9 term. Note 
that for this geometry also in the IR (“confined”) region one has dS/dr ~ 


4 Conclusions 

In this work we have studied the 9 dependence in Witten’s large N c Yang-Mills model using 
the holographic correspondence. 

The model allows to extract the exact 9 dependence of a class of interesting observables, 
like the vacuum energy density, the string tension, the mass of the baryon vertex, the ’t Hooft 
loop, the confinement-deconfinement critical temperature, the mass of the 0 ++ glueball and 
the entanglement entropy. The model shows a common trend of all the mass scales of the 
theory: they get reduced by 9. More precisely, they scale as inverse powers of (1 + 0 2 ) with 

19 We have some different numerical factors w.r.t. fl2l . 

20 Again, we have some different factors w.r.t. [15]. 

21 Modulo a factor of \ due to the two branches of the disconnected surface in the slab case. 
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0 ~ A 4 (0 + 2iik)/N c . Moreover, the model explicitly realizes the expected periodic structure 
of the (T,6) phase diagram, with triple points where first-order transition lines meet, see 
figure [2j It also provides an explicit realization of oblique confinement, with the expected 
relations among the string tensions. 

In the classical gravity regime, valid for strong ’t Hooft coupling (A 4 ^ 1) and large 
N CJ Witten’s model consists of an SU(N C ) gauge theory coupled with adjoint Kaluza-Klein 
massive matter fields. Despite this feature, we have shown that the model shares with pure 
Yang-Mills not only all the relevant IR features, like confinement and the occurrence of a 
mass gap, but also the behavior with 6 of a class of relevant observables. The comparison 
in this case has been done looking at the results obtained on the lattice for small 6 . The 
qualitative agreement we have found in the sign of the first 6 2 corrections let us hope that 
the holographic results, which are exact in 9, can provide useful benchmarks and stimulate 
further lattice analysis on the subleading corrections. 
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A The #-backreacted solution from the Id effective ac¬ 
tion 


Let us look for a regular solution to the action (2.2). 
following ansatz for the lOd metric 


As in m !T8] . let us consider the 


dsl 0 = e 2X dx^dx^ + e 2A dxj 


+ l~s e 


- 2 *dp 2 


+ K e ~ 


'dnl, 


(A.l) 


where p = 0,... 3 label the 4d Minkowski coordinates of the spacetime where the dual 
Yang-Mills held theory is defined, A, A ,ip, v are functions of the radial coordinate p and the 
^-coordinate is compactihed on a circle of length /3 4 = 2tiM^ 1 k . The function ip is related 
to the dilaton </> through the defining equation 

ip = 2 (p — 4A — A — Av . (A.2) 


The transverse part of the metric includes a four-sphere S 4 . 
quantized and proportional to N c , see eq. (2.6). 


The flux of F 4 through S 4 is 
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We will consider the following ansatz for F 2 


F -2 = (IC\ = h dx 4 A dp , 


(A.3) 


where h is a function of the radial variable p and the dot stands for derivative with respect 
to p. 

Implementation of the ansatz above gives the following Id action 

S = V J dp —4A 2 - A 2 - 4z> 2 + ip 2 - I e 4A-A+^+4^2 + v ^ 

V = I2e~ 2u - 2 * - Q 2 c e 4X+ ~ x - 4u ^, (A.4) 

which has to be supported by the zero-energy constraint 

- 4A 2 - A 2 - 4z> 2 + ip 2 - I e 4A-A+ v +4^2 _ v = o _ (A.5) 

Above we have defined (using R 3 = ng s N c l 3 ) 

ft= yk? = Tr- (A ’ 6) 

as the constant arising from the quantized F 4 flux through A 4 and 

1 97r 

(A - 7) 

where V \, 3 is the infinite 4d space-time volume and 4^4 = 87 t 2 /3 is the volume of S 4 . 


A.l Equations of motion and general solutions 

The equation of motion for the field h implies that 

h = - qee - 4X+ ' x -^- 4u = -qee 2X ~ 2 * , 


(A.8) 


where qe is a constant which we will relate to the Yang-Mills 6 angle. The equations of 
motion for the remaining fields, which we re-arrange so that the dilaton (j) appears instead 
of (p, read 

V _ Qc 8A+2A-2 cj> _ © 2A-2 <j> 

' 2 4 

i _ Qc 8A+2A-2 <t> _ _Fe_ 2A-20 
2 4 

$ - .^ e 8A+2A-20 = 

_j_ Qc ^8A+2A—2(/> ^g8A+2A—40+6z^ Qo ^2X—2cf) 
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Interestingly, defining 


7 = 8A + 2A —20, p = 8A + 2A-40 + 6i/, * = 2(A - 0), (A.10) 

we get the following decoupled differential equations 


7 = 4 Qc e 7 , P = 18 e p , x = -2 q 2 e e x , 


(A.ll) 


which can be easily solved in terms of elementary functions 


7 = 


X = 

These in turn imply that 


2 log[ai - e~“ 2P ] - a 2 p + log 
p = -2 log[a 3 - e~ a4p ] - a A p + log 
2 log[a 5 + e~ a&p ] - a 6 p + log 


a\a\ 

Wc 

a 3 al 


a 5 aj 

L % 2 J 


(A- 12 ) 


1 p 

v ~ o (7 + x) + ;: + a 7 + °8 P ■ 
8 6 


(A-13) 

We have 8 integration constants a\,...,as to be hxed imposing IR regularity and the con¬ 
straint (A.5) to begin with. 


Another requirement is that when qg = 0 the solutions read 


1 


Mp) = -4 lo g [1 - e + 4 lo g ^ , 

Mp) = -7 lo g [1 - e~ 3ap ] - lap + j log ^ , 


where we have introduced the parameter 


4 0 L ~ ' J 2 

Mp) = tog [1 - e _3ap ] + \ tog + lo § 9s , 
Mp) = tog [1 - e _3ap ] + 7 log ^ + log y , 


4 
uo 
R 

Uo 

4 0 R 


l. 


(A. 14) 


a = 


MQ c u 0 


Ur 


3R 3 g s l 3 s g 2 s ’ 


(A.15) 


where Uq is the minimal value of the radial variable u used in the main body of this paper. 
As in the qg = 0 case, the variable p is related to the radial variable u by 


3 —3 ap _ -I u o 
1 . .3 > 


u° 


(A.16) 


so that p —> 0 (resp. p —» 00) when u 00 (resp. u —> u 0 ). 
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A.2 Fixing some of the integration constants 

From the expressions above we can deduce the following constraints. First of all, we look 
for a regular confining solution in the 1R. Requiring that e 2A goes to a constant in the 1R 
amounts to choose 

®6 = a 2 • (A.17) 

If in turn we require that, as in the 9 = 0 case, e 2A —y oo in the UV we need to require 


ai = 1. 


(A.18) 


Moreover, requiring e 2u to go to a constant in the IR and to diverge in the UV we find 


1 1 

a% — — -ci2 + -oq , 

4 o 


a 3 = 1 . 


(A.19) 


Requiring e 2A to go to zero in the IR (so that the rr 4 circle shrinks as in the 9 
find the condition 

O4 <C 2fl2• 


0 case) we 
(A. 20 ) 


The zero-energy constraint (A.5) in turn fixes (once all the above conditions are satisfied) 


O4 — (32 • 


(A.21) 


Once all the conditions above are fulfilled the (aq,r) subspace is regular in the deep IR if 


3 2 / 3 a 4 / 3 e 4a T 4 

——— -= -Ml 

V2^l 2 s Q c qe 9 


(A. 22 ) 


A.3 A particular solution 

Since the equations for 7 and p are the same in form as in the qg = 0 case, a first simple 
choice for the yet-to-be fixed integration constants reads 


(32 = 3(3 , 


so that 7 = 70 and p = po. Using the above choice, we get 



Notice that with these choices 


e x(p) 


a 2 a 5 e 3ap 
q 2 e ( 1 + a 5 e 3ap ) 2 


(A.23) 


(A.24) 


(A.25) 





Since this is identified with the running gauge coupling (see e.g. [53]) , and we expect 9 not 
to modify it in the deep UV, we require that 


e x(°) _ e xo(o) _ J_ 
9 2 s 


which sets 


“ 5 - -1 + 9 |f |1 + ’ /1 


9 a 2 g 


2 | ’ 
s 


(A.26) 


(A.27) 


where we have selected the value such that at qg = 0, x — Xo- When qg —>• 0 we have, in fact 

9 a 2 ^ 2 


05 —> 


2 5 

5* 


(A.28) 


so that x —)■ Xo consistently. 


limit, and if we assume <27 to be independent from qg, we can also fix 


If, from (A.22), we want to get the usual relation between w 0 and Mkk in the qg > 0 

(A.29) 


07 = y log g 8 
6 


As for the RR potential, we get 


Ci = hdx 4 = — 
Qe 


05 


a 5 + e~ 3ap 


+ 1 


dx A , 


(A.30) 


where (in order to get a regular solution) we have imposed it to vanish in the deep IR 
(p —y 00) where the the x 4 circle shrinks. The Yang-Mills 6 angle is identified (mod 2n) with 
the integral of Ci along the S X4 circle at p —> 0 


9 + 2nk= - C l = 


2tt 


3a 


27T 2 qg 


S JS; 


X 4 


Mkk Qe ^(1 + 05 ) M KK 3a l s g 2 s ^ A 4gf^ 


(A.31) 


In the expression above fceZ. When qg ~ 0 we get (setting k = 0) 


Qe 

3 al s g 2 ’ 

(A.32) 

1 -A) 

Up ' 

(A.33) 


and so 

W ^ lye" 3 ”'’ = 

Pa Pa 

which is (re-written in terms of the original radial variable u) the known expression one gets 
(for small 6 ) at first order in 6 |3]. 
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Notice that taking the opposite limit qg —> oo is not allowed as we must have 


M < 2 a ^ 


(A.34) 


When qg —> ±(3/2 )ag s we get (setting k = 0) 9 —> ±f3±/(l s g s ) = ±47 t 2 W/A 4 , where A 4 = 


A-K 2 g s N c l s / /3 4 is the ’t Hooft coupling at the KK scale M KK . n On the other hand, the 
validity of the gravity approximation requires 9 <C An 2 N c /X 4 ^ 4 [5J. Thus, while 6 can be of 
order N c as expected, 9/N c must be small. 

A.4 The particular solution in the w-variable 

Let us now see how the particular solution found above reads in terms of the radial coordinate 


u defined, as in the qg = 0 case, in (|A.16|). Let us also define 

so that Iwl < 1. 


2 qe 
w = - 


3 ag s : 


Using (A.31) we get 
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w = 


where, for any k 


0 = 


1±© 2 ’ 

A 4 9 ± 27 rk 
4tt 2 Nr 


(A.35) 

(A.36) 
(A.37) 


Notice that the bound (A.34) is automatically satisfied by any O. 


Using these redefinitions, the particular solution found above turns out to precisely coin¬ 
cide with that written in the main body of the paper and already found in @ 115 ]. 


References 

[1] E. Witten, “Anti-de Sitter space, thermal phase transition, and conEnement in gauge 
theories,” Adv. Theor. Math. Phys. 2 , 505 (1998) |hep-th/980313l]]. 

[2] E. Vicari and H. Panagopoulos, “Theta dependence of SU(N) gauge theories in the 
presence of a topological term,” Phys. Rept. 470 , 93 (2009) |arXiv:0803.1593 [hep-th]]. 

[3] E. Witten, “Theta dependence in the large N limit of four-dimensional gauge theories,” 
Phys. Rev. Lett. 81 , 2862 (1998) |hep-th/9807109|. 

[4] J. L. F. Barbon and A. Pasquinucci, “Aspects of instanton dynamics in AdS / CFT 
duality,” Phys. Lett. B 458 , 288 (1999) |hep-th/99041901. 

22 Notice that when qg — > ±(3/2)og s , a 5 —> 1. 


30 











[5] S. Dubovsky, A. Lawrence and M. M. Roberts, “Axion monodromy in a model of holo¬ 
graphic gluodynamics,” JHEP 1202, 053 (2012) |arXiv:1105.3740 [hep-th]]. 

[6] L. Del Debbio, G. M. Manca, H. Panagoponlos, A. Skouroupathis and E. Vicari, “Theta- 
dependence of the spectrum of SU(N) gauge theories,” JHEP 0606, 005 (2006) |hep- 
th/060304l]j. 

[7] M. D’Elia and F. Negro, lL 9 dependence of the deconfinement temperature in Yang-Mills 
theories,” Phys. Rev. Lett. 109, 072001 (2012) |arXiv:1205.0538 [hep-lat]]. M. D’Elia 
and F. Negro, “Phase diagram of Yang-Mills theories in the presence of a 6 term,” Phys. 
Rev. D 88, no. 3, 034503 (2013) |arXiv: 1306.2919 [hep-lat]]. 

[8] H. Panagopoulos and E. Vicari, “The 4D SU(3) gauge theory with an imaginary 6 
term,” JHEP 1111, 119 (2011) |[arXiv: 1109.6815 [hep-lat]]. 

[9] C. Bonati, M. D’Elia, H. Panagopoulos and E. Vicari, “Change of 6 Dependence in 4D 
SU(N) Gauge Theories Across the Deconfinement Transition,” Phys. Rev. Lett. 110, 
no. 25, 252003 (2013) |arXiv:1301.7640 [hep-lat]]. 

[10] W. Cai, C. Wu and Z. Xiao, “Baryons in the Sakai-Sugimoto model in the D0-D4 back¬ 
ground,” Phys. Rev. D 90, no. 10, 106001 (2014) |arXiv:1410.5549 [hep-th]]. S. w. Li 
and T. Jia, “Matrix model and Holographic Baryons in the D0-D4 background,” 
arXiv: 1506.00068 [hep-th], 

[11] I. R. Klebanov, D. Kutasov and A. Murugan, “Entanglement as a probe of confinement,” 
Nucl. Phys. B 796, 274 (2008) |arXiv:0709.2140 [hep-th]]. 

[12] A. Pakman and A. Parnachev, “Topological Entanglement Entropy and Holography,” 
JHEP 0807, 097 (2008) |arXiv:0805.1891 [hep-th]]. 

[13] C. Wu, Z. Xiao and D. Zhou, “Sakai-Sugimoto model in D0-D4 background,” Phys. 
Rev. D 88, no. 2, 026016 (2013) |arXiv:1304.2111 [hep-th]]. 

[14] D. J. Gross and H. Ooguri, “Aspects of large N gauge theory dynamics as seen by string 
theory,” Phys. Rev. D 58, 106002 (1998) |hep-th/9805129|. 

[15] T. Sakai and S. Sugimoto, “Low energy hadron physics in holographic QCD,” Prog. 
Tlieor. Phys. 113, 843 (2005) jhep-th/0412141]. 

[16] E. Witten, “Large N Chiral Dynamics,” Annals Phys. 128, 363 (1980). 

[17] M. Unsal, “Theta dependence, sign problems and topological interference,” Phys. Rev. 
D 86, 105012 (2012) |arXiv: 1201.6426 [hep-th]]. E. Poppitz, T. Schafer and M. Unsal, 
“Universal mechanism of (semi-classical) deconfinement and theta-dependence for all 


31 




simple groups,” JHEP 1303, 087 (2013) |arXiv:1212.1238 [hep-th]]. M. M. Anber, “0 
dependence of the deconfining phase transition in pure SU(N C ) Yang-Mills theories,” 
Phys. Rev. D 88, no. 8, 085003 (2013) |arXiv:1302.2641 [hep-th]]. 

[18] F. Bigazzi and A. L. Cotrone, “Holographic QCD with Dynamical Flavors,” JHEP 
1501, 104 (2015) [arXiv: 1410.2443 [hep-th]]. 

[19] D. Mateos, R. C. Myers and R. M. Thomson, “Thermodynamics of the brane,” JHEP 
0705, 067 (2007) |hep-th/0701132|. 

[20] G. T. Horowitz and R. C. Myers, “The AdS/CFT Correspondence and a New Posi¬ 
tive Energy Conjecture for General Relativity,” Phy. Rev. D59 026005 (1998), [hep- 
th/9808079]. 

[21] C. Bonati, M. D’Elia, H. Panagopoulos and E. Vicari, “0 dependence of 4D SU(N) 
gauge theories at finite temperature,” PoS LATTICE 2013, 136 (2014) |arXiv:1309.6059 
[hep-lat]]. 

[22] J. M. Maldacena, “Wilson loops in large N held theories,” Phys. Rev. Lett. 80, 4859 
(1998) |hep-th/9803002 . S. J. Rey and J. T. Yee, “Macroscopic strings as heavy quarks 
in large N gauge theory and anti-de Sitter super gravity,” Eur. Phys. J. C 22, 379 (2001) 
|hep-th/98030011. 

[23] F. Bigazzi, A. L. Cotrone, L. Martucci and L. A. Pando Zayas, “Wilson loop, Regge 
trajectory and hadron masses in a Yang-Mills theory from semiclassical strings,” Phys. 
Rev. D 71, 066002 (2005) [hep-th/0409205]. 

[24] A. Brandhuber, N. Itzhaki, J. Sonnenschein and S. Yankielowicz, “Wilson loops, con¬ 
finement, and phase transitions in large N gauge theories from supergravity,” JHEP 
9806, 001 (1998) |hep-th/9803263[. 

[25] C. Schmidhuber, “D-brane actions,” Nucl. Phys. B 467, 146 (1996) |hep-th/96010031. 

[26] C. Csaki, H. Ooguri, Y. Oz and J. Terning, “Glueball mass spectrum from supergravity,” 
JHEP 9901, 017 (1999) [hep-th/980602l(. 

[27] N. R. Constable and R. C. Myers, “Spin two glueballs, positive energy theorems and 
the AdS / CFT correspondence,” JHEP 9910, 037 (1999) jhep-th/9908175]. 

[28] R. C. Brower, S. D. Mathur and C. I. Tan, “Glueball spectrum for QCD from AdS 
supergravity duality,” Nucl. Phys. B 587, 249 (2000) |hep-th/0003115j. 

[29] D. Elander, A. F. Faedo, C. Hoyos, D. Mateos and M. Piai, “Multiscale confining 
dynamics from holographic RG hows,” JHEP 1405, 003 (2014) |arXiv:1312.7160 [hep- 
th]]. 


32 





[30] F. Brunner, D. Parganlija and A. Rebhan, “Glueball Decay Rates in the Witten-Sakai- 
Sugimoto Model,” Phys. Rev. D 91, 106002 (2015) |arXiv:1501.07906 [hep-ph]]. 

[31] F. Brnner and A. Rebhan, “Nonchiral enhancement of scalar glueball decay in the 
Witten-Sakai-Sugimoto model,” arXiv: 1504.05815 [hep-ph], 

[32] G. Gabadadze and A. Iglesias, “On theta dependence of glueballs from AdS / CFT,” 
Phys. Lett. B 609, 167 (2005) |hep-th/041127|. 

[33] E. Witten, “Baryons and branes in anti-de Sitter space,” JHEP 9807, 006 (1998) [hep- 
th/9805112], 

[34] O. Aharony, J. Sonnenschein and S. Yankielowicz, “A Holographic model of decon¬ 
finement and chiral symmetry restoration,” Annals Phys. 322, 1420 (2007) [hep- 
th /0604161]]. 

[35] G. Mandal and T. Morita, “Gregory-Laflamme as the confinement/deconfinement tran¬ 
sition in holographic QCD,” JHEP 1109, 073 (2011) |arXiv:1107.4048 [hep-th]]. 

[36] M. Hanada, Y. Matsuo and T. Morita, “Instanton dynamics in finite temperature QCD 
via holography,” arXiv: 1505.04498 [hep-th], 

[37] O. Bergman and G. Lifschytz, “Holographic U(1)(A) and String Creation,” JHEP 0704, 
043 (2007) |hep-th/0612289|. 

[38] A. Parnachev and A. R. Zhitnitsky, “Phase Transitions, theta Behavior and Instantons 
in QCD and its Holographic Model,” Phys. Rev. D 78, 125002 (2008) |arXiv:0806.1736 
[hep-ph]]. A. S. Gorsky, V. I. Zakharov and A. R. Zhitnitsky, “On Classification of QCD 
defects via holography,” Phys. Rev. D 79, 106003 (2009) jarXiv:0902.1842 [hep-ph]]. 

[39] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from 
AdS/CFT,” Phys. Rev. Lett. 96, 181602 (2006) |hep-th/060300l], 

[40] S. Ryu and T. Takayanagi, “Aspects of Holographic Entanglement Entropy,” JHEP 
0608, 045 (2006) |hep-th/0605073|. 

[41] T. Nishioka, S. Ryu and T. Takayanagi, “Holographic Entanglement Entropy: An 
Overview,” J. Phys. A 42, 504008 (2009) |arXiv:0905.0932 [hep-th]]. 

[42] F. Bigazzi, A. L. Cotrone, C. Nunez and A. Paredes, “Heavy quark potential with 
dynamical flavors: A First order transition,” Phys. Rev. D 78, 114012 (2008) 
|arXiv:0806.1741 [hep-th]]. 


33 






